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• Local emissions of micro-cogeneration system were evaluated. 

• Influence of climatic conditions and control logic was investigated. 

• The micro-cogeneration system was compared with a conventional system. 

• The micro-cogeneration system reduces the NO* emissions. 

• The micro-cogeneration system increases the CO emissions. 
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The increasing diffusion of micro-cogeneration systems is raising the need for studying their environ¬ 
mental impact in order to assess their sustainability. The adoption of the systems for the combined 
production of heat and power may provide a significant reduction of global impact in terms of carbon 
dioxide emissions with respect to the separate production of electricity and heat. However, a compre¬ 
hensive environmental evaluation of this technology should take into account as well the impact due to 
the presence of plants spread over the territory that could increase the local pollution, in particular due 
to nitrogen oxides and carbon monoxide, and thus could worsen the local air quality. 

In this paper the nitrogen oxides and carbon monoxide emissions of a residential building-integrated 
micro-cogeneration system were evaluated; a 6.0 kW e i natural gas fuelled internal combustion engine- 
based micro-cogeneration unit was coupled with a multi-family house compliant with the transmittance 
values suggested by the Italian Law. The analyses were carried out by using the whole building simu¬ 
lation software TRNSYS upon varying the city where the building is located (four Italian cities repre¬ 
sentative of different climatic regions were considered) as well as the control logic of cogeneration device 
(electric and thermal load-following strategies). The simulated performance of the proposed system was 
compared with those of a conventional system composed of a natural gas-fired boiler (for thermal en¬ 
ergy production) and a power plant mix connected to the central electric grid (for electricity production) 
in order to assess the suitability of the cogeneration-based system in reducing the local emissions. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Micro-cogeneration (MCHP) is the combined production of heat 
and power with electric output lower than 50 kW e i [1 . The heat 
produced is used for space and water heating and possibly for 
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cooling, while the electricity is used within the building or fed into 
the grid [2,3]. The “core” of MCHP systems is a prime mover based 
on different technologies: reciprocating internal combustion en¬ 
gine (ICE), reciprocating external combustion Stirling engine (SE), 
micro-turbine (MT), and fuel cell (FC); among the others technol¬ 
ogies, reciprocating internal combustion engine-based micro¬ 
cogeneration systems are usually the best choice for small scale 
cogeneration applications thanks to their robust nature and reli¬ 
ability, high thermal efficiency, fuel flexibility, reduced pollutant 
emissions, low noise and vibration, small installation space [2,3 . 
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In the last decade, technological innovations and a changing 
economic and regulatory environment have resulted in a renewed 
interest for distributed generation [4]. At present, environmental 
policies or concerns are probably the major driving force for the 
demand of distributed generation in Europe with a number of 
research programmes around the world, aimed at improving the 
performance of small-scale energy systems, as well as reducing the 
related emissions [5 . 

Several studies demonstrated the effectiveness of MCHP sys¬ 
tems in reducing the primary energy consumption with respect to 
the conventional separate production of heat and electricity with a 
percentage difference varying from 10% to 30%, depending on the 
size and efficiency of the cogeneration units [5-7 . The reduction in 
the fuel consumption from such MCHP systems brings a corre¬ 
sponding reduction of global emissions of CO 2 in a first approxi¬ 
mation similar to the amount of energy saving 8-10]. However, 
whilst basically all CHP plants provide primary energy saving with 
respect to the heat-and-electricity separate production, and thus 
potential CO 2 emission reduction at a global level, the scenario may 
be not so good for what concerns the local emissions 11,12 . In fact 
the diffusion of distributed cogeneration within urban areas, where 
air quality standards are quite stringent, brings environmental 
concerns on a local level due to the emission of hazardous air 
pollutants such as NO x , CO, SO x , particulate matter, unburned hy¬ 
drocarbons, and so on. Indeed, in urban contexts dispersion in the 
atmosphere of pollutants from small-scale generators sited among 
buildings may be more difficult than, for instance, for large-scale 
power plants with high stacks. In addition, also due to the high 
population density, there is a number of relative weak receptors 
(elderly and sick people, children, etc.), with other potential im¬ 
pacts of pollutant emissions referred to ecosystems, monuments, 
and so forth. A further critical point is represented by the already 
high background emission level mostly due to road traffic pollution. 
As a consequence, air quality standards and emission limits can be 
quite stringent in urban areas, and the local impact of MCHP sys¬ 
tems over the environment is an issue to analyse with care; 
therefore a comprehensive evaluation of micro-cogeneration sys¬ 
tems should not consider the global environmental impact due to 
the CO 2 emissions only, but it should also take into account the local 
environmental impact associated to the presence of plants spread 
over the territory that could increase the local pollution, and thus 
could worsen the local air quality. 

The most relevant hazardous pollutants emitted by natural gas- 
fired equipments are CO and NO x (often subject to binding envi¬ 
ronmental constraints [13 ). NO x formation is the one giving most 
concerns from a regulatory point of view, being its toxic effects 
turning up for concentrations 10-20 times lower than for CO 13]. 
Subsequently, NO x reduction represents also a research field in 
which several efforts have been put over the last years [14 . 

At the moment few works are available in literature concerning 
the evaluation of the local impact associated to cogeneration sys¬ 
tems. Canova et al. [11 ] characterized both a 60 l<W e i micro-turbine 
and a 180 l<W e i internal combustion engine fed with natural gas in 
terms of air pollutants (in particular, NO x , CO and non-methane 
organic compounds). The results from the case study showed that 
in general the micro-turbine always performs better than the in¬ 
ternal combustion engine in terms of NO x (but when it is operated 
at partial load the CO levels increase dramatically), while the CO 
emissions from the ICE are only negligibly affected by the engine 
load. Mancarella and Chicco 12] presented a general methodology 
to assess the emission impact due to adoption of MCHP systems in 
urban areas under general operating conditions (including, in 
particular, partial-load operation). They applied this methodology 
by considering two different MTs available on the market (with 30 
and 60 l<W e i as nominal electric output, respectively) fuelled with 


natural gas. Their analyses highlighted how the full-load emission 
characterization can be unsuitable to estimate the environmental 
impact due to actual operating conditions; critical situations have 
emerged in the specific case for what concerns high NO x emissions 
at low load (for the 30 l<W e i MT) and high CO emissions mainly at 
medium-low loading condition. Pehnt [15 investigated the po¬ 
tential environmental effects of a deployment of micro-cogenera¬ 
tion on different levels. He noticed that whereas for fuel cells and 
Stirling engines emissions of air pollutants are extremely low, 
reciprocating engines emit more significant amounts of NO x , CO, 
and hydrocarbons. Furthermore, the emission factors of recipro¬ 
cating engines depend heavily on operation characteristics, age and 
maintenance of the systems. In addition to investigating the 
emissions side, analysis of the concentration of pollutants in a 
residential area supplied by reciprocating engines was carried out: 
the analysis showed that for selected, rather unfavourable, mete¬ 
orological conditions, with a flat topography and urban housing 
structure, the additional level of NO x due to the engines does not 
create severe additional environmental impacts. In a previous work 

[16] the authors estimated the NO x and CO emissions of a resi¬ 
dential building-integrated micro-cogeneration system during the 
winter. A commercially available 6.0 kW e i internal combustion 
engine based-MCHP unit fuelled with natural gas [17] was used to 
satisfy both thermal and electric requirements of a multi-family 
house composed of three floors compliant with the transmittance 
values of both walls and windows suggested by the Italian Law 18]. 
The whole system was analysed by means of the dynamic whole¬ 
building simulation software TRNSYS [19 (1 min simulation time 
step was used) in the case of the building located in the city of 
Naples (south of Italy) with the MCHP device operating under 
thermal load-following logic. In the current work the evaluation 
was extended in a wider range of operating conditions. In partic¬ 
ular, in order to evaluate the influence of climatic conditions on the 
system environmental performance, a comprehensive sensitivity 
analysis was performed by considering the building located into 
four Italian cities representative of different climatic regions: 

> Palermo (latitude: 38° 07' North; longitude: 13° 20' East); 

> Napoli (latitude: 40° 50' North; longitude: 14° 15' East); 

> Roma (latitude: 41° 53' North; longitude: 12° 29' East); 

> Milano (latitude: 45° 28' North; longitude: 9° 10' East). 

In addition to the thermal load-following logic, the operation of 
the MCHP unit under electric load-control strategy was also 
investigated with the main aim to assess the influence of the MCHP 
control logic on the local environmental impact. 

The emissions of NO x and CO calculated for the proposed system 
were compared with those of a conventional system composed of a 
natural gas-fired boiler (for thermal energy production) and a po¬ 
wer plant mix connected to the national central grid (for electricity 
production) in order to assess the suitability of the MCHP unit 
based-system in comparison to the one based on separate energy 
production from an environmental point of view. 

The analyses performed in Refs. 11,12,15] were performed by 
considering data on local emissions provided by the manufacturers, 
while in this paper the local emissions of the MCHP unit were 
calculated on the basis of the data measured at both the Advanced 
Engine Technology Ltd. in Ottawa (Canada) and the Canadian 
Centre for Housing Technology twin research homes 20] during 
full and part-load operation of the above-mentioned MCHP device 

[17] . Taking into consideration that often little attention is paid to 
the emission worsening that could be brought by consistent oper¬ 
ation of micro-cogeneration devices at partial loads (this could lead 
to biased environmental assessment of the systems evaluated only 
on the basis of the full-load performance, whereas load-tracking 
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operation can frequently occur), in this paper the part-load char¬ 
acteristics of the cogeneration unit was fully considered by using a 
detailed dynamic model [21,22] calibrated and validated by the 
authors based on experimental results [23,24 . The local emissions 
associated to both the boiler and central electric grid were calcu¬ 
lated by using the standard emission factors approach suggested in 
Ref. [12]. 

In Sections 2 and 3 the proposed system is described in detail, 
while Section 4 provides the main characteristics of the system 
assumed as reference. The applied methodology and the main re¬ 
sults obtained from the simulations are reported and commented 
in the Section 5. 

2. Description of the proposed system 

Fig. 1 shows the scheme of the system configuration analysed in 
this study. The whole plant basically consists of a micro¬ 
cogeneration device (MCHP), an auxiliary boiler (BOILER), a com¬ 
bined storage tank (TANK) with three internal heat exchangers 
(IHE1, IHE2, IHE3), a plate-fin heat exchanger (PHE), three pumps 
(PI, P2, P3), three thermostats (Tl, T2, T3), two 3-way valves (VI, 
V2), two flow diverters (Dl, D2), and a group of fan-coils installed 
within the building. 

Four different main circuits can be easily distinguished in Fig. 1 : 

1) micro-cogeneration device circuit; 

2) natural gas-fired boiler circuit; 

3) space heating circuit; 

4) domestic hot water circuit. 

Both natural gas and electricity lines are also highlighted in 

Fig. 1. 

The system under investigation is devoted to satisfy the space 
heating sensible load, the domestic hot water requirements and the 
electric demand related to the heating season of a multi-family 
house composed of three floors. 

The heating purposes and domestic hot water production are 
satisfied by heating up the water contained within the tank. The 
MCHP system and natural gas-fired boiler co-operate in order to 
guarantee a given water temperature level (55 °C) within the hot 
water storage. The group of fan-coils is supplied by the combined 
tank. According to the European Standard EN12831:2003 [25], 
20 °C was assumed as set-point indoor air temperature. Domestic 
hot water is produced by means of the internal heat exchanger IHE3 
located in the tank; in the case of the temperature of the water 


exiting IHE3 is lower than 45 °C, the required additional heat is 
provided by the auxiliary natural gas-fired boiler and transferred to 
the domestic hot water by means of the plate heat exchanger PHE. 

The electricity generated by the cogeneration device is used 
directly into the building for the lighting systems, domestic appli¬ 
ances (such as vacuum cleaner, dishwasher, washing machine, PC, 
TVs, fridge), auxiliaries (pumps, fans, etc) or is exported to the 
electric grid; the electric grid is also used to cover peak demand. 

3. Components and loads of the proposed system 

Each component of the whole system was simulated by using 
the software TRNSYS [19 , where each physical piece is modelled 
with a component (named “Type”) that is a FORTRAN source code 
model. In the following sections the main features of the Types 
used in this study for simulating each plant component are re¬ 
ported. Additional details can be found in Ref. 26]. 

3.2. Combined tank 

The combined tank for both heating purposes and domestic hot 
water production was modelled in this study by means of the 
Type60f included in TRNSYS library. This model allows to calculate 
thermal stratification in the component: in this study ten temper¬ 
ature levels (nodes) were used in the tank; a uniform tank loss 
coefficient per unit area equal to 3.0 kj/hm 2 I< was assumed [27 . A 
vertical cylindrical hot water storage unit with one flow inlet (IF in 
Fig. 1) and one flow outlet (OF in Fig. 1) was considered: the cold 
water coming from the building enters the tank through the flow 
inlet IF, while hot water going towards the fan-coils installed within 
the building exits the tank through the flow outlet OF. The tank was 
equipped with three internal heat exchangers: the hot water 
coming from the MCHP unit flows through the internal heat 
exchanger located in lower part of the tank (IHE1 in Fig. 1 ); the hot 
water coming from the natural gas-fired boiler goes towards the 
internal heat exchanger located in upper part of the tank (IHE2 in 
Fig. 1 ); the third internal heat exchanger (IHE3 in Fig. 1 ) allows heat 
to be extracted for domestic hot water production. The main 
characteristics of the tank considered in this study are reported in 
Table 1 [28]. 

3.2. MCHP unit 

A reciprocating internal combustion engine-based micro¬ 
cogeneration unit commercialized by AISIN SEIKI company 17] was 


MULTI-FAMILY HOUSE 


BOILER CIRCUIT 
MCHP CIRCUIT 
DHW CIRCUIT 
SPACE HEATING CIRCUIT 
NATURAL GAS LINE 
ELECTRICITY LINE 



Fig. 1 . Scheme of the system configuration analysed in this study. 
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Table 1 

Main characteristics of the analysed hot water storage [28]. 


Tank volume (m 3 ) 0.855 

Tank height (m) 2.040 

Height of flow inlet (m) 0.125 

Height of flow outlet (m) 1.691 

Temperature levels used in the tank 10 

Tank loss coefficient (kj/hm 2 K) 3.0 

Destratification conductivity (kj/hm K) 0 

Number of internal heat exchangers 3 

Height of thermostat controlling the MCHP (m) 0.816 

Height of thermostat controlling the boiler (m) 1.836 

The fraction of glycol contained in the IHE1 by volume (%) 0 

IHE1 inside diameter (m) 0.027 

IHE1 outside diameter (m) 0.029 

IHE1 total surface area (m 2 ) 3.1 

IHE1 length (m) 37.0 

IHE1 wall conductivity (W/m K) 45 

IHE1 material conductivity (W/m K) 45 

Height of IHE1 inlet (m) 0.813 

Height of IHE1 outlet (m) 0.144 

The fraction of glycol contained in the IHE2 by volume (%) 0 

IHE2 inside diameter (m) 0.026 

IHE2 outside diameter (m) 0.028 

IHE2 total surface area (m 2 ) 2.5 

IHE2 length (m) 30.5 

IHE2 wall conductivity (W/m K) 45 

IHE2 material conductivity (W/m K) 45 

Height of IHE2 inlet (m) 1.521 

Height of IHE2 outlet (m) 1.048 

The fraction of glycol contained in the IHE3 by volume (%) 0 

IHE3 inside diameter (m) 0.023 

IHE3 outside diameter (m) 0.025 

IHE3 total surface area (m 2 ) 7.8 

IHE3 length (m) 106.5 

IHE3 wall conductivity (W/m K) 16 

IHE3 material conductivity (W/m K) 16 

Height of IHE3 inlet (m) 0.108 

Height of IHE3 outlet (m) 1.703 


investigated in this study. The main characteristics of this device 
are reported in Table 2. 

Thermal power is recovered from the exhaust gases and the 
engine jacket of the micro-cogenerator by means of a water-glycol 
mixture and transferred to the water within the storage tank by 
using water as coolant. The MCHP device under investigation can 
operate under both electric and thermal load-following control 
strategies. Both operation logics were simulated in this study: 

Table 2 

Main characteristics of the MCHP unit [17]. 


AISIN SEIKI unit 


Model 

GECC60A2 

Engine type 

Reciprocating internal combustion engine, 
water cooled, 4 strokes, 3 cylinders 

Displacement (cm 3 ) 

952 

Speed revolution (rpm) 

1600-1800 

Fuel 

Natural gas, LPG 

Generator type 

Permanent — magnet type, synchronous 
generator 16 poles 

Rated electric output (kW) 

6 

Rated thermal output (kW) 

11.7 

Operating sound at 

54 

1.0 m distance and 

1.5 m height (dB) 

Electric efficiency" 1 at 

28.8 

maximum load (%) 

Thermal efficiency at 

56.2 

maximum load (%) 


a Ratio between electric output and primary power input. 
b Ratio between heat supplied to the end-user and primary power input. 


a) in the case of thermal load-following operation, the system 
operates according to the external temperature signal coming 
from the thermostat T1 (Fig. 1) placed on the tank near the 
entrance of the internal heat exchanger IHE1: when this tem¬ 
perature is lower than the set-point value (55 °C), the unit 
provides its maximum electric and thermal outputs; when the 
water temperature in the storage exceeds the set-point value, 
the unit is turned off. The electricity generation is a by-product 
and any unused excess of electricity is sent to the power line, 
that is also used to cover the peak load; 

b) in the case of electric load-following operation, the cogeneration 
system produces electricity by modulating the power supplied 
through the inverter on the basis of the electric demand. The 
MCHP system operates only in the case of the electric load is 
larger than 300 W: if the electric load is lower than 300 W the 
cogeneration unit does not start, while the MCHP system pro¬ 
duces its maximum power in the case of the electric load is 
larger than its maximum electric output. Heat is also available 
during unit operation, even if its level depends on the system 
electric output. 

The AISIN SEIKI unit has an internal control to keep the internal 
water-glycol temperature lower than a given value (about 70 °C), in 
order to prevent conditions endangering the engine block. This 
particular behaviour of the cogeneration device was taken into 
account in this study in the case of electric load-following opera¬ 
tion: in the case of the temperature of the coolant coming out from 
the cogeneration device becomes higher than 70 °C, the recovered 
heat is dissipated into the external air through an air-to-water heat 
exchanger installed within the unit. 

The Italian Law does not limit the daily operation time of 
cogeneration systems: as a consequence, in this study the MCHP 
unit was allowed to operate up to 24 h per day depending on 
electric/thermal requirements. The cogeneration unit was simu¬ 
lated by using the detailed dynamic model ; 21,22] developed 
within the Annex 42 of the International Energy Agency. This 
model was calibrated and validated for the AISIN SEIKI unit [17] on 
the basis of several laboratory tests performed by the authors 
[23,24] at the Built Environment Control Laboratory of the Second 
University of Naples. During the experiments the transient and 
steady-state operation of the cogeneration system was experi¬ 
mentally investigated upon varying the electric level, the coolant 
temperature, as well as the coolant flow rate. The experimental 
tests showed no significant variations of the electric efficiency of 
the micro-cogeneration device between transient and stationary 
regime; on the contrary, the values of thermal efficiency measured 
during transient operation were significantly lower than those 
measured during steady-state operation: this highlighted that 
neglecting the influence of thermal transition may imply a relevant 
overestimation of thermal energy supplied by the unit to the end- 
user, mainly during system operation at low electric levels. The 
analysis of the tests also indicated that both electrical and thermal 
efficiencies during steady-state operation are not significantly 
affected by both coolant temperature and coolant flow rate, while 
they are strongly dependent from the electric output of the unit: in 
particular, the experiments showed a thermal efficiency increasing 
when electric load reduces, and an electric efficiency increasing 
with the electrical power. The measurements were compared with 
the simulation data in order to empirically validate the calibrated 
Annex 42 model; the data measured during a 24 h dynamic test 
performed by applying a realistic daily load profile representing the 
Italian domestic non-HVAC demand profile for a multi-family 
house of five dwellings were used: the predicted fuel use, the 
electric output and heat recovery over the whole duration of the 
data set differed from measured values by -6.1%, -0.2% and -1.0%, 
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respectively. Taking into account the uncertainties associated to 
both measuring instruments and calibration procedure [23,24], it 
can be stated that the results of the performed validation activity 
demonstrated the suitability of the Annex 42 model in predicting 
the performance of the MCHP unit; in fact similar errors were 
detected in Refs. [29,30]. 

3.3. Boiler 

A 20.0 kW t h natural gas-fired boiler [31 was considered in this 
study. The boiler was modelled in TRNSYS by using the Type6 
included in TRNSYS library. The auxiliary heater is activated only in 
the case of (i) the water temperature within the tank (sensed by the 
thermostat T2 in Fig. 1 ) is lower than a given set-point value (55 °C), 
or when (ii) the domestic hot water temperature at the outlet of the 
internal heat exchanger IHE3 (sensed by the thermostat T3 in Fig. 1 ) 
is lower than 45 °C. In the case of temperature measured by ther¬ 
mostat T3 being lower than 45 °C, the full boiler output is utilized 
for hot water through the plate heat exchanger PHE. 

According to the manufacturer data [31 ], 9.5 1/min as water flow 
rate was considered in this study. 

The manufacturer [31 suggests 0.924 as boiler efficiency for a 
part load ratio equal to 30% and 0.927 as boiler efficiency in the case 
of operation at rated output. Therefore in this paper the boiler ef¬ 
ficiency ? 7 B ’ SP was evaluated by considering a linear interpolation of 
the manufacturer data: 

rj B SP = 0.924 + 0.000214- (p B ; sp - 6.0j (1) 

Exhaust and combustion losses rates were not considered as 
negligible in this work: therefore the primary power consumption 
of the boiler was evaluated as follows: 

P p SP = P t B h SP /^ ,SP ( 2 ) 

According to the Italian Law [32], Italy is divided in 6 different 
climatic zones (named A, B, C, D, E, F) based on the Heating Degree 
Days (HDD) index [33]. A maximum daily operation time for the 
heating system depending on the climatic zone is specified in 
Ref. [32]; as a consequence, the boiler was allowed to operate only 
during the time intervals specified in Table 3 as a function of both 
the city and day (week day and weekend day are distinguished). 
Out of these periods the boiler was turned off for heating purposes; 
the boiler is allowed to operate up to 24 h per day if thermal energy 
is required for producing domestic hot water at 45 °C. 

3.4. Building characteristics and loads 

The geometrical layout of the building investigated in this report 
is basically a multiplication of a single family house type building 
geometry. The main geometrical characteristics of each floor and the 
building orientation are shown in Fig. 2. All three floors have the 
same useable floor area (96.0 m 2 ), while the net height of each single 
floor is 3.0 m. Five windows were considered for each single floor. 


Starting from 1st January 2010, the Italian Law [18] specifies the 
threshold values of transmittance for both walls and windows of 
renovated buildings depending on (i) the climatic zone where the 
building is located, and (ii) the wall type (external wall, ground and 
roof). These threshold values of transmittance Uii m for property 
renovations located in Palermo, Napoli, Roma e Milano are specified 
for walls and windows, respectively, in Tables 4and 5. 

In Table 4 the characteristics of building walls considered in this 
study are also specified: in particular, the internal (ft;) and external 
( h e ) convective and radiative heat transfer coefficients of walls, wall 
transmittance (U wa n), layers material, thickness (s), density (p), 
specific heat (c) and thermal conductivity (A) are reported for each 
wall type as well as for each city. As it can be derived from Table 4, 
each wall is compliant to the requirements of Italian Law 18], with 
transmittance values equated to the given threshold values: this 
result was obtained by adjusting the building insulation upon 
varying the climatic zone. The values of hi and h e were derived by 
the European Standard EN ISO 6946 [34], while the values of p, c 
and A were defined based on Italian Standard UNI 10351:1994/EC 
[35]. 

The main characteristics of windows are reported in Table 5: in 
particular, the internal (hi) and external (h e ) convective and radia¬ 
tive heat transfer coefficients of windows, frame transmittance 
tAframe, glazing transmittance U g i az i ng , window type, glazing width, 
spacing width, spacing gas, ratio between frame area Af rame and 
window (glazing + frame) area A W i n d 0 w are specified for each city. 
As it can be derived from Table 5, windows transmittance values are 
compliant to the requirements of Italian Law [18 with trans¬ 
mittance values equated to the given threshold values; values of 
fAframe and U W i n dow were calculated based on the European Standard 
EN ISO 10077-1 [36]; all windows are provided with shadings, 
whose zone and wall reflection coefficients were considered equal 
to 0.5. According to the Italian Law [37 , the area of each window 
was defined as the 12.5% of the ground area of the room where the 
window is located. 

The Italian Law [18] specifies the duration of the heating season 
depending on the Italian region where the plant is located: ac¬ 
cording to the limits specified in Ref. [18 , in this study the 
following periods of simulations were considered: 

i^ Palermo: from December 1st to March 31st; 

i^ Napoli: from November 15th to March 31st; 

i^ Roma: from November 1st to April 15th; 

i^ Milano: from October 15th to April 15th. 

A maximum daily operation time for the heating system 
depending on the climatic zone is specified in Ref. 18]. As a 
consequence, the indoor air temperature was controlled only dur¬ 
ing the time intervals specified in Table 3: out of these time in¬ 
tervals the indoor air temperature was not controlled. 

Heat coming from occupants, personal computers and lighting 
systems was assumed to contribute to the internal gains of the 
building. Fig. 3 shows the number of occupants of each single floor 
and the occupants-related sensible heat gain of each single floor as 


Table 3 

Time intervals during which the boiler can operate for heating purposes. 

Week day Weekend day 

Palermo Napoli Roma Milano Palermo Napoli (climatic zone C, HDD = 1034) Roma Milano 

(climatic zone B, (climatic zone C, (climatic zone D, (climatic zone E, (climatic zone B, (climatic zone D, (climatic zone E, 

HDD = 751) HDD = 1034) HDD = 1415) HDD = 2404) HDD = 751) HDD = 1415) HDD = 2404) 

7:30-10:00 7:30-11:00 6:30-11:00 6:30-12:00 8:30-10:00 8:30-11:00 8:30-12:00 8:30-13:30 

16:00-21:30 15:00-21:30 15:00-22:30 14:30-23:00 16:00-22:30 15:00-22:30 14:00-22:30 14:30-23:30 
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12.0 m 

Windows area: 

W1 = 2.23 nr 
W2 = 1.04 nr 
W3 = 2.75 nr 
W4 = 1.82 nr 
W5 = 3.00 nr 


Fig. 2. Main geometrical characteristics 
house. 



s 


a single flat of the whole multi-family 


a function of the time during a week day (Fig. 3a) and a weekend 
day (Fig. 3b). Sensible heat coming from each occupant was 
assumed equal to 75.0 W: this values is suggested by the Standard 
ISO 7730 [38] in the case of light work/typing as degree of activity. 

Fig. 4 shows the sensible heat flux related to the lighting ap¬ 
pliances installed within the rooms of each single flat. Light sources 
with an installed total electric capacity of 294.0 W were consid¬ 
ered; thermal power coming from each lighting system was 
assumed equal to the 75% of its nominal electric capacity. Sensible 


heat produced by lighting systems was considered as transferred to 
the indoor air by both radiation (70%) and convection (30%) [39 . 

Concerning the infiltration, the European Standard 
EN12831:2003 [25] was used to calculate the air exchange rate, i.e. 
the number of interior volume air changes that occur per hour, 
induced by wind and stack effect on the building envelope: based 
on this European Standard, the air exchange rate was assumed 
equal to 0.28 h -1 . 

The internal gains coming from occupants, personal computers 
and lighting appliances for the entire building were generated 
through the superposition of three of the single flat profiles. No 
difference among the four cities was considered in terms of internal 
gains. 

A group of fan-coils, supplied by the hot water storage, was 
installed into the building in order to balance the space heating 
sensible load. A temperature difference of 5 °C was assumed for the 
hot water crossing the fan-coils, with a variable-flow pump con¬ 
trolling the fan-coils outlet temperature. 

The Type56a (included in TRNSYS library) was used for model¬ 
ling the building envelope, indoor air set-point temperature, infil¬ 
tration and internal gains. Typel5-6 (included in TRNSYS library) 
was used for reading the external EnergyPlus weather data files 
[40] of each city. A specific EnergyPlus weather data file was 
considered for each city. 

Fig. 5 reports the whole building thermal load-duration diagram 
with the heat-demand values sorted in descending order for each 
city under consideration. As it can be derived from Fig. 5, space 
heating demand has a quite different duration depending on the 
city under investigation: it is equal to about 940 h for Palermo, 
1370 h for Napoli, 1860 h for Roma, and 2550 h for Milano. In 
addition it can be noticed that the building sensible thermal load is 
lower than the nominal thermal output of the MCHP unit (11.7 kW) 
for the most part of the system operation time whatever the city 
under consideration is. The maximum thermal load is always 
31.7 kW: this value was set in TRNSYS simulations as the sum of the 


Table 4 

Characteristics of building walls. 


Wall type 

u lim (W/m 2 I<) 

hi (W/m 2 K) 

h e (W/m 2 K) 

Uwaii (W/m 2 K) 

Layer materials 

s(m) 

P (kg/m 3 ) 

c (kJ/Kg k) 

A (W/m K) 

External wall 

Palermo = 0.48 

7.7 

25.0 

Palermo = 0.48 

Gypsum plaster 

0.02 

900 

1.00 

0.250 


Napoli = 0.40 



Napoli = 0.40 

Bricks 

0.33 

1800 

0.84 

0.720 


Roma = 0.36 



Roma = 0.36 

Expanded Polyester 

Palermo = 0.053 

20 

1.20 

0.041 


Milano = 0.34 



Milano = 0.34 

Gypsum plaster 

Napoli = 0.070 
Roma = 0.082 

Milano = 0.088 
0.02 

900 

1.00 

0.250 

Roof 

Palermo = 0.38 

10.0 

25.0 

Palermo = 0.38 

Gypsum plaster 

0.015 

1800 

1.00 

0.700 


Napoli = 0.38 



Napoli = 0.38 

Cement and bricks 

0.300 

1000 

0.84 

0.300 


Roma = 0.32 



Roma = 0.32 

Concrete 

0.084 

1900 

0.84 

1.060 


Milano = 0.30 



Milano = 0.30 

Bitumen 

0.010 

1200 

1.00 

0.170 






Expanded polyester 

Palermo = 0.052 
Napoli = 0.052 
Roma = 0.072 
Milano = 0.080 

20 

1.20 

0.041 






Sand 

0.050 

1700 

0.84 

1.400 






Lime and cement adhesive 

0.010 

1500 

1.00 

0.720 






Tiles 

0.010 

2300 

0.84 

0.720 

Ground 

Palermo = 0.49 

5.9 

25.0 

Palermo = 0.49 

Tiles 

0.010 

2300 

0.84 

0.720 


Napoli = 0.42 



Napoli = 0.42 

Lime and cement adhesive 

0.010 

1500 

1.00 

0.720 


Roma = 0.36 



Roma = 0.36 

Sand 

0.050 

1700 

0.84 

1.400 


Milano = 0.33 



Milano = 0.33 

Expanded polyester 

Palermo = 0.026 
Napoli = 0.040 
Roma = 0.057 

Milano = 0.067 

20 

1.20 

0.041 






Bitumen 

0.005 

1200 

1.00 

0.170 






Concrete 

0.077 

1900 

0.84 

1.060 






Bricks 

0.300 

1000 

0.84 

0.300 






Gypsum plaster 

0.015 

1800 

1.00 

0.700 
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Table 5 

Characteristics of windows. 



tflim 

(W/m 2 K) 

hi 

(W/m 2 K) 

h e 

(W/m 2 K) 

tfframe 

(W/m 2 K) 

tfglazing 

(W/m 2 K) 

Window type 

Glazing 
width (mm) 

Spacing 
width (mm) 

Spacing 

gas 

^frame/^window (^) 

Palermo 

3.0 

7.7 

25.0 

3.96 

2.83 

Double glazing 

4 

16 

Argon 

15 

Napoli 

2.6 



2.94 

2.54 


6 




Roma 

2.4 



7.90 

1.43 


4 

15 



Milano 

2.2 



6.56 

1.53 
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Fig. 3. Occupants-related sensible heat gains for a single flat during week day (a) and weekend day (b) whatever the city is. 


thermal output of the boiler of the proposed system (20.0 kW) and 
the maximum thermal output of the MCHP unit (11.7 kW). 

Fig. 6 shows the values of thermal energy required for heating 
purposes of the whole building per unit area as a function of both 
the city and the simulation period. According to this figure, it is 
worth noting that, whatever the period under investigation is, the 
maximum and the minimum amounts of thermal energy are 
required in the case of the building located in Milano and Palermo, 
respectively; the thermal energy required by the building located in 
Napoli is always larger than that one used by the building located in 
Roma. 

In this study the domestic hot water demand profile with an 
average basic load of 200 1/day in the time scale of 1 min specified 


within IEA-SHC Task 26 [41,42 was used for estimating the domestic 
hot water demand of each single flat. Fig. 7a shows the daily domestic 
hot water demand of each single flat for three typical days (December 
31 st, January 31 st and March 15th), while Fig. b reports the hot water 
demand-duration diagram of each single flat as a function of the city, 
with the domestic hot water flow rate values sorted in descending 
order. In this paper the domestic hot water demand profile of the 
entire building was generated through the superposition of three of 
the single flat profiles. The temperature of cold water coming from 
the mains was derived by considering the distribution profile defined 
by the Typel5-6 (included inTRNSYS library). 

Fig. 8a highlights the daily electric demand profile resulting 
from the operation of both lighting systems and other domestic 



Time (hh:mm) 


Fig. 4. Daily heat flux profile associated to the lighting appliances of a single floor whatever the city is. 
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Fig. 5. Thermal load-duration diagram with the heat-demand values sorted in descending order for each city under consideration. 
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Fig. 6. Thermal energy required for building heating purposes per unit area as a function of the city and the simulation period. 




Fig. 7. Single flat daily domestic hot water demand profile during 3 typical days (a) and single flat domestic hot water duration curve (b). 
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appliances (such as vacuum cleaner, dishwasher, washing machine, 
PC, TVs, fridge) assumed in this study for each single flat. The 
electric consumptions of each appliance were derived from 
Ref. [43 . The daily electric demand profile of the whole building 
was defined through the superposition of three of the single flat 
profiles reported in Fig. 8a. The electric demand profile considered 
in this study corresponds to an electric consumption for the whole 
building equal to about 109.7 Wh/m 2 per day. Fig. 8b shows the 
electric load-duration diagram (with the electric-demand values 
sorted in descending order) for each city under consideration. This 
figure shows that the building requires an electric power larger 
than the nominal electric output of the micro-cogeneration unit 
(6.0 kW) for a time interval ranging from 52 h (Palermo) and 79 h 
(Milano). 

The electric consumption of the pumps and the fan of heat 
dissipator were also taken into account during the simulations. 

4. Description of the reference system 

The main focus of this study is comparing the environmental 
performance of the proposed plant (SP) based on a micro¬ 
cogeneration device with those of a conventional system (CS) 
based on separate energy production. In the following a conven¬ 
tional system composed of a natural gas-fired boiler (for thermal 
energy production) and a power plant mix connected to the electric 
grid (for electricity production) is considered. The comparison was 
performed by assuming for the reference system the same electric 
and thermal energy output of the proposed system. 

A 32.0 kW t h natural gas-fired boiler [31] was selected for the 
thermal energy production within the reference system. The boiler 
manufacturer [31 suggests a value of 0.911 as boiler efficiency for a 
part load ratio equal to 30% and a value of 0.937 as boiler efficiency 
in the case of operation at rated output; in this study the boiler 
efficiency of the conventional system n] B ^ cs was evaluated by 
considering a linear interpolation of the manufacturer data: 

rj B CS = 0.911 + 0.001067 - (p B ,; cs - 10.44) (3) 

The exhaust and combustion losses rates were not considered as 
negligible in this study. Therefore the primary power consumption 
of the boiler considered for the conventional system was evaluated 
as follows: 


P BCS = P t B [; cs /4 B ’ CS (4) 

Concerning the efficiency of the power plant mix connected to 
the national electric grid ri PPl a figure of 0.461 was assumed: this 
value represents the power plant average efficiency in Italy, 
including transmission losses [44 . 

5. Methodology, results and discussion 

In this section NO* and CO emissions of the proposed system 
were evaluated and compared with those associated to the con¬ 
ventional system assumed as reference. 

Regardless of the excess air, any combustion process is charac¬ 
terized by a certain amount of by-products such as, in particular, CO 
not oxidized into CO 2 . In general, the emissions of this pollutant are 
not easy to calculate, since they depend on several specific com¬ 
bustion kinetic characteristics, such as temperature, oxygen excess, 
mixing aspects, and so forth. In addition, the emission levels are 
commonly kept at a minimum at the design stage, and their in¬ 
crease during operation is an indicator of lower combustion effi¬ 
ciency (for instance, due to bad fuel-air mixing or bad 
maintenance) [11,13,45]. The NO x formation mainly depends on 
temperature in the combustion chamber, retention time in the 
primary combustion zone, as well as the relative air-to-fuel ratio 
(i.e. the ratio of the actual air to fuel ratio to the stoichiometric one) 
that affects not only the NO* production, but also the CO emissions 
[11 . NO x and CO emissions of the micro-cogeneration system [18] 
under investigation in this paper were not measured by the au¬ 
thors, but they were derived from the experimental data measured 
at both the Advanced Engine Technology Ltd. and the Canadian 
Centre for Housing Technology twin research homes [20]. The 
measurement results were obtained by using a Horiba PG250 
portable gas analyzer [46]; this analyzer was connected to the 
MCHP unit by means of a sampling probe inserted approximately 
50 cm into the flue outlet. Emissions of NO x and CO were tested by 
operating the cogeneration device under electric load-following 
control strategy, with the electric load varied in stepwise manner 
from about 2.5 kW to about 6 kW; levels of NO x and CO emissions 
were recorded at a frequency of 1 s -1 . Testing was performed in two 
phases: in the first phase cold start testing was performed, while 
the second phase was carried out for considering the case of warm- 
start. The cold weather tests began with a cold start, where the 
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Fig. 8. Daily electric demand profile for each single flat (a) and building electric load-duration diagrams (b). 
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engine had been kept outside for 2.5 h in -17 °C ambient air prior 
to start-up; cold start testing was done with the unit installed in an 
environmental chamber for soaking in temperatures down 
to -40 °C. The warm weather testing phase was conducted with the 
cogeneration unit installed outdoors. In this paper only the values 
measured during the second phase were considered. In addition it 
must be highlighted that the unit experimentally investigated in 
Ref. [20] was the model commercialized by the AISIN SEIKI com¬ 
pany 17] within the non-European Union market: this model is not 
equipped with the catalytic converter for abating the CO emissions; 
the model commercialized by the AISIN SEIKI company [17] for the 
European Union market is equipped with the catalytic converter, 
and therefore it is characterized by lower emissions of CO. 

Fig. 9 shows the emissions of both NO* and CO as a function of 
the electric output of the MCHP system measured during steady- 
state operation in warm weather conditions as reported in 
Ref. [20 . As it can be derived from this figure, steady-state NO x 
levels were about 60 ppm at 6 kW and 18 ppm at 2.5 kW; CO 
emissions of about 560 ppm at full power and 800 ppm at 2.5 kW 
were observed. The NO x and CO emissions were strongly correlated 
with electric load within the range examined. Emissions were 
measured by Lombardi et al. [20] also during start-up, shutdown 
and load transients; the highest CO emissions were observed dur¬ 
ing start-up, where levels of 1400 ppm were obtained briefly 
approximately 5 min after engine firing; the NO x levels did not 
noticeably spike during start-up. In the following only the emis¬ 
sions measured during steady-state operation (reported in Fig. 9) 
were taken into account. 

In this study the emissions of both NO x and CO due to the MCHP 
operation were evaluated by using the following formulae: 


investigated in Ref. [20 . Trend lines expressed by Eqs. (5) and (6) 
are also reported in Fig. 9. 

The emissions of both NO x and CO due to the MCHP operation 
expressed in kilograms were obtained starting from Eqs. (5) and (6) 
by means of the following formulae: 

= eK , ™ P /10 6 -M NOl /M exh -m fuer (l + «)• STS (7) 

m MC H p = e MCHP /10 6 -M co /M exh ■ m fuel • (1 + a) • STS (8) 

where M is the molar mass of the chemical compound 
(M NOv = 44 kg/kmol, Mqo = 28 kg/kmol), M ex h is the average molar 
mass of exhaust gases (assumed equal to 29 kg/kmol), rh fuel is the 
fuel mass flow rate consumed by the MCHP unit (calculated by 
assuming the lower heating value of natural gas equal to 49,599 kj/ 
kg [47]), STS is the Simulation Time Step used in this study (equal to 
1 min). The parameter a represents the air to fuel ratio; its assumed 
value (17.12 kg a i r /kgf ue i) was obtained by considering the following 
complete stoichiometric methane combustion: 

CH 4 + 20 2 + 7.52N 2 = C0 2 + 2H 2 0 + 7.52N 2 (9) 

The emissions of NO x and CO due to both the operation of the 
natural gas-fired boiler of the proposed system and the production 
of electricity purchased from the central electric grid were calcu¬ 
lated by using the standard emission factors approach suggested in 
Ref. 12]. According to this approach, the mass of a given pollutant x 
emitted while producing the energy output E can be worked out as: 

m x = u p -E (10) 


eft 1 ™ 13 = 0.5346-(P e ^ CHP ) 3 -2.8313-(P e ^ CHP ) 2 

+ 5.5675-P^ chp +13.707 (5) 

e MCHP = _ 4 779 . (pMCHP ) 3 + 66.843- (P“ CHP ) 2 

-361.02-P^ chp +1362.5 (6) 

where P ( )j ,CHP is the electric output of the cogeneration device. The 
values provided by Eqs. (5) and (6) are expressed in parts per 
million. Both Eqs. (5) and (6) were defined by correlating the 
experimental data reported in Fig. 9 with a third order polynomial 
function; a coefficient of determination equal to one was obtained 
in both cases within the range of electric power experimentally 


where u\ is the energy output-based emission factor that is the 
specific emissions of x per unit of E. As a consequence, the mass of 
both NO x and CO emitted by the natural gas-fired boiler of the 
proposed system were calculated as follows: 



(ii) 


m co P = (12) 

E E 

where u N p 0 and u c p 0 are, respectively, the NO x and CO emission 
factors associated to the natural gas consumption, E^ SP is the 
thermal energy produced by the boiler of the proposed system, and 
? 7 b ’ SP is the efficiency of the boiler of the proposed system. The 
values of ? 7 B ’ SP were calculated by means of Eq. (1). 
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Fig. 9. Steady-state emissions during warm weather testing 20]. 
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In this study the NO* and CO emission factors associated to the 
natural gas consumption were assumed equal to 0.2 gNO x /l<Wh p 
and 0.04 gCO/kWh p , respectively, according to the values suggested 
in Ref. [12 : these values represent the equivalent NO* and CO 
emissions for a primary energy unit consumed by burning out 
natural gas. 

The mass of both NO x and CO emitted for producing the elec¬ 
tricity purchased from the central electric grid were obtained by 
using the following formulae: 


m 


grid 

NO, 


Eel cgrid 
u NO, c e\ 


(13) 


<o d = u co' £ d id (14) 

The average NO x and CO emissions for producing electricity 
depend upon the technology mix adopted. Thus, the emissions are 
significantly lower where the generation system is based upon 
renewable sources such as hydroelectric power. Italy has a power 
generation system with no nuclear power and with about 20% of 
hydro-electricity; the main part is thermal-based, with a mix of 
fuels including coal, oil and natural gas. According to the values 
suggested in Ref. [12 , in this paper the NO x and CO emission factors 
for electricity production were assumed equal to 0.5 gNO x /kWh e i 
and 0.3 gCO/kWh e i, respectively: these value represents the 
equivalent NO x and CO emissions for an electric energy unit 
consumed according to the Italian scenario. 

Table 6 highlights the MCHP operation time, the electric energy 
£MC Hp produced by the cogeneration device, the thermal energy 
£M CHp SU ppii e d to the water within the tank by the MCHP unit, the 
thermal energy £^ sp produced by the boiler of the proposed sys¬ 
tem, and the electric energy Ef { ld purchased from the central grid. 
These parameters are reported as a function of both the city under 
consideration and the MCHP unit control logic. This table shows 
that, in the case of thermal load-following logic, the thermal energy 
supplied by the MCHP system to the tank is 4.7 -=-7.7 times larger 
than those produced by the boiler, while the electric output of the 
AISIN SEIKI unit is 1.02.9 times greater than the electric energy 
purchased from the central national grid; in the case of electric 
load-control strategy, the values of £{^ chp are 0.42.6 times the 

values of £^ SP , while E ^ cHP is 11.7-^13.1 times greater than £^ nd . 

Table 7 shows the mass of both NO x and CO emitted by each 
single components of the proposed system (the MCHP device, the 
natural gas-fired boiler as well as the electric grid) during the 
whole heating season as a function of both the city under consid¬ 
eration and the MCHP unit control logic. This table highlights that: 

> under electric load following logic, the emissions of NO x asso¬ 
ciated to the MCHP unit are about 1.12.1 times larger than 
those associated to the boiler, and about 6.6 -=-7.3 times greater 
than those related to the electric grid; 

> under electric load following logic, the emissions of CO associ¬ 
ated to the MCHP unit are about 107.9-^197.9 times larger than 
those associated to the boiler and about 326.5^-356.1 times 
greater than those related to the electric grid; 

> under thermal load following logic, the emissions of NO x asso¬ 
ciated to the MCHP unit are about 3.9-f-5.6 times larger than 
those associated to the boiler, and about 0.82.2 times the 
values related to the electric grid; 

> under thermal load following logic, the emissions of CO asso¬ 
ciated to the MCHP unit are about 168.9-^-261.0 times larger 
than those associated to the boiler and about 19.1 -^44.6 times 
greater than those related to the electric grid; 


Table 6 

Energy flows entering and exiting the proposed system during the whole heating 
season. 




Electric load¬ 
following operation 

Thermal load¬ 
following operation 

Palermo 

MCHP operation 
time (h) 

2136 

1180 


£ e T HP (kJ) 

2.02E + 07 

1.38E + 07 


EK CHP (kJ) 

C P (kJ) 

3.75E + 07 

3.79E + 07 


1.51E + 07 

4.90E + 06 


£?r (kj) 

1.72E + 06 

1.34E + 07 

Napoli 

MCHP operation 
time (h) 

2541 

1639 


E“ CHP (kJ) 

2.47E + 07 

2.44E + 07 


E“ CHP (kj) 

4 SP (kJ) 

4.33E + 07 

5.68E + 07 


3.09E + 07 

1.20E + 07 


(kJ) 

1.88E + 06 

1.25E + 07 

Roma 

MCHP operation 
time (h) 

3132 

1932 


E e T HP (kJ) 

3.02E + 07 

2.21E + 07 


r“ CHP (kj) 

4 SP (kJ) 

5.38E + 07 

6.01 E + 07 


2.08E + 07 

8.1 IE + 06 


£fr (kj) 

2.37E + 06 

1.93E + 07 

Milano 

MCHP operation 
time (h) 

3563 

2822 


E“ CHP (kJ) 

3.57E + 07 

4.28E + 07 


E“ CHP (kj) 

C P (kJ) 

1.66E + 07 

9.55E + 07 


4.35E + 07 

1.76E + 07 


Ef, nd (kJ) 

2.86E + 06 

1.49E + 07 


> in comparison to the electric load-control strategy, the thermal 
load-following operation is characterized by 1.31.4 times 
larger total NO x emissions and 0.50.7 times lower total CO 
levels. 

As above-mentioned, the emissions of NO x and CO associated to 
the proposed system (SP) were compared with those of the con¬ 
ventional system (CS) described in the Section 4. The comparison 
was performed by means of the following indicators: 



(15) 


Table 7 

Mass of both NO, and CO emitted by each single components of the proposed system 
during the whole heating season. 


Electric load- Thermal load¬ 
following operation following operation 


Palermo 


Napoli 


Roma 


Milano 


MCHP 


mzA"' (kg) 


m 


NO. 

grid 


(kg) 


m No:( k g) 
m m hP (kg) 

m co, P (kg) 

<o d (kg) 

nC“P( k g) 

m No!(kg) 

<a(kg) 

m m hP (kg) 

<o P (kg) 

<d (kg) 

mKK HP (kg) 
m No! (kg) 
<a(kg) 
m m hP (kg) 

mgr (kg) 

grid, 

CO 

jjjMCHP rVnr') 

BSf, (Kg) 


(kg) 


m 


NO. 

grid 


(kg) 


< ( k §) 
m [$p P (k § ) 


m 

m 


CO. 

grid 

CO 


(kg) 

(kg) 


1.58E + 00 
8.36E-01 
2.39E - 01 
3.12E + 01 
1.67E - 01 
9.54E - 02 
1.91E + 00 
1.72E + 00 
2.61E - 01 
3.71E + 01 
3.44E - 01 
1.04E - 01 
2.37E + 00 
1.16E + 00 
3.30E - 01 
4.57E + 01 
2.31E - 01 
1.32E - 01 
2.77E + 00 
2.42E + 00 
3.98E - 01 
5.22E + 01 
4.83E - 01 
1.59E - 01 


1.52E + 00 
2.72E-01 
1.86E + 00 
1.42E + 01 
5.45E - 02 
7.46E - 01 
2.59E + 00 
6.64E - 01 
1.73E + 00 
2.24E + 01 
1.33E - 01 
6.92E - 01 
2.39E + 00 
4.51E - 01 
2.68E + 00 
2.14E + 01 
9.02E - 02 
1.07E + 00 
4.46E + 00 
9.80E - 01 
2.07E + 00 
3.69E + 01 
1.96E - 01 
8.28E - 01 
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Rco = m co/ m co 


,cs 


(16) 


where the N0 X and CO emissions of both the proposed system and 
the conventional system are calculated as follows: 


m NO x — m NO x 


MCHP , m B,SP grid 
+ " l NO x m NO x 


(17) 


,SP 


m co = m co 


MCHP , m B,SP grid 
T m CO " l co 


(18) 


in 


cs 

NO 


X 



( cMCHP 
(^th 




+ U 


Bel 

NOv 


£^ chp + £®[ id ) 


(19) 


111 


CS 

CO 



( cMCHP 




+ U 


E e i 

CO 


^MCHP 



( 20 ) 


A value of R NOx and R C o lower than one means that the pro¬ 
posed system allows to reduce the emissions of NO* and CO, 
respectively, with respect to the system based on separate energy 
production. 

Figs. 10 and 11 report the values of R NOa and R c o as a function of 
both the city under consideration and the MCHP unit control logic 
under thermal and electric load-control strategy, respectively. 
These figures highlight that: 


whatever the city and MCHP control logic are, the proposed 
system allows for reducing the NO x emissions in comparison 
to the conventional system: the NO x levels associated to the 
proposed system are 0.44-^0.59 times lower than those 
related to the reference system; 
i^ whatever the city and MCHP control logic are, the CO emis¬ 
sions associated to the conventional system are substantially 
negligible with respect to those of the proposed system: the 
CO levels associated to the proposed system are 7.37 17.47 
times larger than those related to the reference system. 
About this point, it should be highlighted that this results is 
mainly due to the fact that in this study the MCHP unit 
commercialized by the AISIN SEIKI company 17 without the 
catalytic converter was considered; 
u* whatever the city, the reduction of the NO x emissions is more 
pronounced in the case of electric load-following logic in 
comparison to the thermal load-control strategy; 



Fig. 10. R NOx and R C o under thermal load-following operation as a function of the city. 



Fig. 11. R NOa and R C o under electric load-following operation as a function of the city. 


is* whatever the city, the thermal load-following logic allows to 
obtain lower levels of CO with respect to the electric load- 
control strategy; 

u* whatever the MCHP control logic is, the influence of climatic 
conditions is quite relevant in terms of both NO x and CO 
emissions: 

a) under thermal load-following operation: the value of R NOa 
ranges from 0.59 when the building is located in Palermo 
up to 0.52 if Milano is the city under consideration; the 
value of Rco is 7.37 when the building is located in Roma 
and 8.50 in the case of Milano; 

b) under electric load-following logic: locating the building 
in Roma allows to achieve the minimum value of R NOa 
( 0.44), while the maximum value (0.50) of R NOa is ob¬ 
tained when the building is located in Napoli or Milano; 
the values of R c o range from 16.05 (Milano) up to 17.47 
(Roma). 

The results showed that, if compared to the conventional 
system assumed as reference, the proposed system is able to 
reduce the NO x emissions, but it is characterized by larger CO 
emissions; therefore the simulation data highlighted that the 
MCHP unit is not completely suitable in reducing the local 
emissions. 

However it can be noticed that, according to the experimental 
measurements, the CO emissions of the micro-cogeneration device 
decrease significantly with the electric load; therefore the thermal 
load-following control strategy should be considered instead of the 
electric-load following logic in order to mitigate the CO emissions 
associated to the proposed system. In addition it should be 
underlined that the MCHP unit analysed in this study is not 
equipped with the catalytic converter for abating the CO emissions 
because the model commercialized by the AISIN SEIKI company 
[17] within the non-European Union market was investigated; the 
model commercialized by the AISIN SEIKI company [17] for the 
European Union market is equipped with the catalytic converter 
and it is characterized by significantly lower emissions of CO. 
Therefore in the future the emissions of carbon monoxide associ¬ 
ated to the device commercialized by the AISIN SEIKI company [17] 
for the European Union market should be also measured in detail 
and compared to the CO emissions associated to the conventional 
systems in order to obtain a comprehensive indication about the 
suitability of the MCHP device from the point of view of local 
emissions. 
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6. Conclusions 

In this paper the NO x and CO emissions of a residential building- 
integrated micro-cogeneration system were evaluated and 
compared with those of a conventional system composed of a natural 
gas-fired boiler (for thermal energy production) and a power plant 
mix connected to the central electric grid (for electricity production). A 
6.0 kW e i reciprocating internal combustion engine fuelled with nat¬ 
ural gas 17] was considered; its NO* and CO emissions were derived 
from experimental measurements carried out by Lombardi et al. 20]. 

The analyses were carried out by using the whole building 
simulation software TRNSYS upon varying the city where the res¬ 
idential building is located as well as the MCHP control logic. 

The simulation results showed that: 

> in comparison to electric load-control strategy, thermal load- 
control strategy is characterized by lower CO emissions, but it 
is less favourable in terms of NO* levels; 

> whatever the city and MCHP control logic are, the proposed 
system allows for reducing the NO* emissions in comparison to 
the conventional system, while the CO emissions associated to 
the conventional system are substantially negligible with 
respect to those of the proposed system; 

> whatever the city, the reduction of the NO x emissions is more 
pronounced in the case of electric load-following logic in com¬ 
parison to the thermal load-control strategy; 

> whatever the MCHP control logic is, the influence of climatic 
conditions is quite relevant in terms of both NO x and CO emissions. 

In the future the CO emissions associated to the device 
commercialized by the AISIN SEIKI company [17] for the European 
Union market (equipped with the catalytic converter) should be 
also evaluated and compared to the carbon monoxide emissions 
associated to the conventional systems based on separate energy 
production in order to get a final indication about the suitability of 
the MCHP device from the point of view of local emissions. 

Of course the numerical results presented in the case study are 
representative of the specific situation, but, even if the specific 
results are not meant to be generalized, the synthetic and simple 
methodology followed in this study is highly suitable for running 
effective analyses with indicative or conventionally established 
emission data, regardless of the specific application. This aspect 
could be particularly useful to formulate emission policies for 
distributed energy systems. 
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Nomenclature 


Latin letters 

A: area (m 2 ) 

c: specific heat (kj/kg K) 

e: emission (ppm) 

E: energy (kWh) 

h: heat transfer coefficient (W/m 2 K) 
M: molar mass (kg/kmol) 
m; mass (kg) 
m; mass flow rate (kg/s) 

P: power (kW) 
s: thickness (m) 

STS: simulation time step (s) 

U ; transmittance (W/m 2 K) 
u: emission factor (g/kWh) 

Greeks 

a: air to fuel ratio (kg air /kgf U ei) 
p: efficiency (-) 


A: thermal conductivity (W/m K) 
p: density (kg/m 3 ) 

Acronyms 

CO: carbon monoxide 
FC: fuel cell 

HDD: Heating Degree Days 

ICE: reciprocating internal combustion engine 

IF: flow inlet 

IHE: internal heat exchanger 

MCHP: micro-combined heat and power 

MT: micro-turbine 

OF; flow outlet 

NO x : nitrogen oxides 

R: ratio between emissions of proposed system and conventional system 
SE: reciprocating external combustion Stirling engine 

Superscripts 
B: boiler 

CS: conventional system 
E: energy 

E p : primary energy 
E e i'■ electric energy 

grid: related to the central electric grid 
MCHP: micro-combined heat and power 
SP: proposed system 

Subscripts 
air: air 

CO; carbon monoxide 
e: external 
el: electric 
exh: exhaust gases 

frame: related to the frame of the window 
fuel: natural gas 

glazing: related to the glass of the window 

i: internal 

lim: limit value 

NO x : nitrogen oxides 

p: primary 

PP: power plant 

th: thermal 

window: related to the whole window 
x: generic pollutant 


